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ABSTRACT 

We present an analysis of ~390 ksec of data of the Z source GX 340+0 taken during 24 obser- 
vations with the Rossi X -ray Timing Explorer satellite. We report the discovery of a new broad 
component in the power spectra. The frequency of this component varied between 9 and 14 
Hz, and remained close to half that of the horizontal branch quasi-periodic oscillations (HBO). 
Its rms amplitude was consistent with being constant around ~5%, while its FWHM increased 
with frequency from 7 to 18 Hz. If this sub-HBO component is the fundamental frequency, 
then the HBO and its second harmonic are the second and fourth harmonic component, while 
the third harmonic was not detected. This is similar to what was recently found for the black 
hole candidate XTE J1550-564. The profiles of both the horizontal- and the normal branch 
quasi-periodic oscillation peaks were asymmetric when they were strongest. We describe this in 
terms of a shoulder component at the high frequency side of the quasi-periodic oscillation peak, 
whose rms amplitudes were approximately constant at ~4% and ~3%, respectively. The peak 
separation between the twin kHz quasi-periodic oscillations was consistent with being constant 
at 339+8 Hz but a trend similar to that seen in, e.g. Sco X-l could not be excluded. We discuss 
our results within the framework of the various models which have been proposed for the kHz 
QPOs and low frequency peaks. 



Subject headings: accretion, accretion disks - 
— X-rays: stars 

1. INTRODUCTION 

GX 340+0 is a bright low-mass X-ray binary (LMXB) 
and a Z source (Hasinger & van der Klis 1989). 
The Z-shaped track traced out by Z sources in the 
X-ray color-color diagram or hardness-intensity dia- 
gram (HID) is divided into three branches: the hor- 
izontal branch (HB), the normal branch (NB), and 
the flaring branch (FB). The power spectral proper- 
ties and the HID of GX 340+0 were previously de- 
scribed by van Paradijs et al. (1988) and Kuulk- 
ers &: van der Klis (1996) using data obtained with 
the EXOSAT satellite, by Penninx et al. (1991) us- 
ing data obtained with the Ginga satellite, and by 
Jonker et al. (1998) using data obtained with the 
Rossi X -ray Timing Explorer (RXTE) satellite. An 
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extra branch trailing the FB in the HID has been 
described by Penninx et al. (1991) and Jonker et al. 
(1998). When the source is on the HB or on the upper 
part of the NB, quasi-periodic oscillations (QPOs) 
occur with frequencies varying from 20-50 Hz: the 
horizontal branch quasi-periodic oscillations or HBOs 
(Penninx et al. 1991; Kuulkers & van der Klis 1996; 
Jonker et al. 1998). Second harmonics of these HBOs 
were detected by Kuulkers & van der Klis (1996) and 
Jonker et al. (1998) in the frequency range 73-76 Hz 
and 38-69 Hz, respectively. In the middle of the NB, 
van Paradijs et al. (1988) found normal branch oscil- 
lations (NBOs) with a frequency of 5.6 Hz. Recently, 
Jonker et al. (1998) discovered twin kHz QPOs in 
GX 340+0. These QPOs have now been seen in all 
six originally identified Z sources (Sco X-l, van der 

1 Astronomical Institute "Anton Pannekoek", University of Amsterdam, and Center for High-Energy Astrophysics, Kruislaan 403, 
1098 SJ Amsterdam; peterj@astro.uva.nl, michiel@astro.uva.nl, homan@astro.uva.nl, mariano@astro.uva.nl, ecford@astro.uva.nl 
2 MIT, Center for Space Research, Cambridge, MA 02139, Chandra Fellow; rudy@space.mit.edu 
3 University of Alabama, Huntsville 

4 Facultad de Ciencias Astronomicas y Geofisicas, Universidad Nacional de La Plata, Paseo del Bosque S/N, 1900 La Plata, 
Argentina 

5 Space Research Organization Netherlands, Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands; e.kuulkers@sron.nl 
"Astronomical Institute, Utrecht University, P.O. Box 80000, 3507 TA Utrecht, The Netherlands 

7 Department of Physics and Astronomy, University of Illinois at Urbana-Champaign, Urbana, IL 61801; f-lamb@uiuc.edu 



2 



Klis et al. 1996; Cyg X-2, Wijnands et al 1998a; 
GX 17+2, Wijnands et al. 1997b; GX 349+2, Zhang 
et al. 1998; GX 340+0, Jonker et al. 1998; GX 5-1, 
Wijnands et al. 1998b; see van der Klis 1997, 1999 
for reviews), but not in Cir X— 1, which combines Z 
source and atoll source characteristics (Oosterbroek 
et al. 1995; Shirey, Bradt, and Levine 1999; see also 
Psaltis, Belloni, & van der Klis 1999). 
In the other class of LMXBs, the atoll sources 
(Hasinger & van der Klis 1989), kHz QPOs are ob- 
served as well (see van der Klis 1997, 1999 for re- 
views). Recently, also HBO-like features have been 
identified in a number of atoll sources (4U 1728- 
34, Strohmayer et al. 1996, Ford & van der Klis 
1998, Di Salvo et al. 1999; GX13+1, Homan et al. 
1998; 4U 1735-44, Wijnands et al. 1998c; 4U 1705- 
44, Ford, van der Klis, & Kaaret 1998; 4U 1915-05, 
Boirin et al. 1999; 4U 0614+09, van Straaten et al. 
1999; see Psaltis, Belloni, & van der Klis 1999 for 
a summary). Furthermore, at the highest inferred 
mass accretion rates, QPOs with frequencies near 6 
Hz have been discovered in the atoll sources 4U 1820- 
30 (Wijnands, van der Klis, & Rijkhorst 1999c), and 
XTE J1806-246 (Wijnands & van der Klis 1998e, 
1999b; Revnivtsev, Borozdin, & Emelyanov 1999), 
which might have a similar origin as the Z source 
NBOs. 

At low mass accretion rates the power spectra of black 
hole candidates, atoll, and Z sources show similar 
characteristics (van der Klis 1994a, b). Wijnands & 
van der Klis (1999a) found that the break frequency 
of the broken power law which describes the broad- 
band power spectrum, correlates well with the fre- 
quency of peaked noise components (and sometimes 
narrow QPO peaks) observed in atoll sources (includ- 
ing the millisecond X-ray pulsar SAX J1808. 4-3658; 
Wijnands & van der Klis 1998d, Chakrabarty & Mor- 
gan 1998), and black hole candidates. The Z sources 
followed a slightly different correlation. In a similar 
analysis, Psaltis, Belloni, & van der Klis (1999) have 
pointed out correlations between the frequencies of 
some of these QPOs and other noise components in 
atoll sources, Z sources, and black hole candidates, 
which suggests these phenomena may be closely re- 
lated across these various source types, or at least de- 
pend on a third phenomenon in the same manner. Be- 
cause of these correlations, models describing the kHz 
QPOs which also predict QPOs or noise components 
in the low-frequency part of the power spectrum can 
be tested by investigating this low-frequency part. 
In this paper, we study the full power spectral range 
of the bright LMXB and Z source GX 340+0 in or- 
der to further investigate the similarities between the 
atoll sources and the Z sources, and to help con- 
strain models concerning the formation of the differ- 



ent QPOs. We report on the discovery of two new 
components in the power spectra of GX 340+0 with 
frequencies less than 40 Hz when the source is on the 
left part of the HB. We also discuss the properties of 
the NBO, and those of the kHz QPOs. 

2. OBSERVATIONS AND ANALYSIS 

The Z source GX 340+0 was observed 24 times in 
1997 and 1998 with the proportional counter array 
(PCA; Jahoda et al. 1996) on board the RXTE 
satellite (Bradt, Rothschild & Swank 1993). A log 
of the observations is presented in Table |]. Part of 
these data (observations 1, 9-18) was used by Jonker 
et al. (1998) in the discovery of the kHz QPOs in 
GX 340+0. The total amount of good data obtained 
was ~390 ksec. During ~19% of the time only 3 or 4 
of the 5 PCA detectors were active. 
The data were obtained in various modes, of which 
the Standard 1 and Standard 2 modes were always 
active. The Standard 1 mode has a time resolution 
of 1/8 s in one energy band (2-60 keV). The Stan- 
dard 2 mode has a time resolution of 16 s and the 
effective 2-60 keV PCA energy range is covered by 
129 energy channels. In addition, high time resolu- 
tion data (with a resolution of 244 /j,s or better for 
the 2-5.0 keV band and with a resolution of 122 /zs 
or better for the 5.0-60 keV range) were obtained for 
all observations. 

For all observations except observation 1, which had 
only 4 broad energy bands, and observation 22, for 
which technical problems with the data occurred, we 
computed power spectra in five broad energy bands 
(2-5.0, 5.0-6.4, 6.4-8.6, 8.6-13.0, 13.0-60 keV) with a 
Nyquist frequency of 256 Hz dividing the data in seg- 
ments of 16 s length each. We also computed power 
spectra for all observations using 16 s data segments 
in one combined broad energy band ranging from 5.0- 
60 keV with a Nyquist frequency of 4096 Hz. 
To characterize the properties of the low-frequency 
part (1/16-256 Hz) of the power spectrum we exper- 
imented with several fit functions (see Section ||) 
but finally settled on a fit function that consisted 
of the sum of a constant to represent the Poisson 
noise, one to four Lorentzians describing the QPOs, 
an exponentially cut-off power law component, P cx 
v~ a exp{—u jvcut) to describe the low frequency noise 
(LFN), and a power law component to represent the 
very low frequency noise (VLFN) when the source 
was on the NB. 

To describe the high frequency part (128 to 4096 Hz 
or 256 to 4096 Hz) of the power spectrum we used a 
fit function which consisted of the sum of a constant 
and a broad sinusoid to represent the dead-time mod- 
ified Poisson noise (Zhang et al. 1995), one or two 
Lorentzian peaks to represent the kHz QPOs, and 
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Table 1 
Log of the observations. 



Number 


Observation 


Date Sz 


Total on source 




ID 


Start time (UTC) 


observing time (ksec.) 


1 


20054-04-01-00 


1997-04-17 13:26:21 


19.8 


2 


20059-01-01-00 


1997-06-06 06:05:07 


34.7 




20059-01-01-01 


1Q97-06-06 21-39-06 


8 1 


4 


20059-01-01-02 


1997-06-07 11:15:05 


22.1 


5 


20059-01-01-03 


1997-06-07 23:48:56 


21.6 


6 


20059-01-01-04 


1997-06-08 07:51:04 


22.9 


7 


20059-01-01-05 


1997-06-09 00:09:03 


17.5 


8 


20059-01-01-06 


1997-06-10 01:22:46 


22.0 


9 


20053-05-01-00 


1997-09-21 01:04:06 


17.5 


10 


20053-05-01-01 


1997-09-23 04:09:30 


11.5 


11 


20053-05-01-02 


1997-09-25 01:30:29 


8.4 


12 


20053-05-01-03 


1997-09-25 09:37:51 


19.3 


13 


20053-05-02-00 


1997-11-01 22:38:58 


9.5 


14 


20053-05-02-01 


1997-11-02 03:32:07 


9.0 


15 


20053-05-02-02 


1997-11-02 19:42:00 


12.7 


16 


20053-05-02-03 


1997-11-03 01:50:07 


13.9 


17 


20053-05-02-04 


1997-11-04 01:59:34 


11.1 


18 


20053-05-02-05 


1997-11-04 16:18:27 


7.2 


19 


30040-04-01-00 


1998-11-13 23:52:00 


16.7 


20 


30040-04-01-01 


1998-11-14 13:55:00 


17.3 


21 


30040-04-01-02 


1998-11-14 21:03:00 


28.1 


22 


30040-04-01-03 


1998-11-15 13:48:00 


17.1 


23 


30040-04-01-04 


1998-11-15 20:57:00 


17.0 


24 


30040-04-01-05 


1998-11-15 09:53:00 


2.6 
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sometimes a power law to fit the lowest frequency 
part (< 150 Hz). The PCA setting concerning the 
very large event window (Zhang et al. 1995; van der 
Klis et al. 1997) was set to 55 /is. Therefore, its effect 
on the Poisson noise was small and it could be incor- 
porated into the broad sinusoid. The errors on the 
fit parameters were determined using Ax 2 =1.0 (la 
single parameter). The 95% confidence upper limits 
were determined using Ax 2 =2.71. 
We used the Standard 2 data to compute hardnesses 
and intensities from the three detectors that were al- 
ways active. Figure [l] shows three HIDs; one (A) 
for observations 1 and 9-18 combined (data set A), 
one (B) for observation 2-8 combined (data set B), 
and one (C) for observation 19-24 combined (data 
set C). The observations were subdivided in this way 
because the hard vertex, defined as the HB-NB in- 
tersection, is at higher intensities in data set C than 
in data set A. The hard vertex of data set B falls at 
an intermediate intensity level. 

We assigned a value to each power spectrum accord- 
ing to the position of the source along the Z track us- 
ing the S z parameterization (Dieters & van der Klis 
1999, Wijnands et al. 1997) applied to each HID 
separately. In this parametrization, the hard vertex 
(defined as the HB-NB intersection) is assigned the 
value S z = 1.0 and the soft vertex (defined as the 
NB-FB intersection) is assigned S z = 2.0. Thus, the 
distance between the hard and soft vertex defines the 
length scale along each branch. Since for HID C we 
only observed part of the Z track, we used the posi- 
tion of the soft vertex of HID A in HID C. From the 
fact that the soft vertex of the HID B was consistent 
with that from HID A, we conclude that the error 
introduced by this is small. 

The shifts in the position of the hard vertex prevented 
us from selecting the power spectra according to their 
position in an HID of all data combined. We selected 
the power spectra according to the S z value in each of 
the three separate Z tracks, since Jonker et al. (1998) 
showed that for GX 340+0 the frequency of the HBO 
is better correlated to the position of the source rela- 
tive to the instantaneous Z track than to its position 
in terms of coordinates in the HID. The power spec- 
tra corresponding to each S z interval were averaged. 
However, employing this method yielded artificially 
broadened HBO peaks, and sometimes the HBO pro- 
file even displayed double peaks. The reason for this 
is that in a typical S z selection interval of 0.05 the 
dispersion in HBO frequencies well exceeds the sta- 
tistical one, as shown in Figure 0. While the relation 
between S z and HBO frequency is roughly linear, the 
spread is large. 

For this reason, when the HBO was detectable in the 
5.0-60 keV power spectra, we selected those power 



spectra according to HBO frequency rather than on 
S z value. In practice, this was possible for all data 
on the HB. To determine the energy dependence of 
the components, the 2-5.0, 5.0-6.4, 6.4-8.6, 8.6-13.0, 
13.0-60 keV power spectra were selected according 
to the frequency of the HBO peak in the 2-60 keV 
power spectrum, when detectable. 
The HBO frequency selection proceeded as follows. 
For each observation we constructed a dynamical 
power spectrum using the 5-60 keV or 2-60 keV data 
(see above), showing the time evolution of the power 
spectra (see Fig. ||). Using this method, we were able 
to trace the HBO frequency in each observation as 
a function of time. We determined the maximum 
power in 0.5 Hz bins over a range of 2 Hz around the 
manually identified QPO frequency for each power 
spectrum, and adopted the frequency at which this 
maximum occurred as the HBO frequency in that 
power spectrum. This was done for each observation 
in which the HBO could be detected. 
The 18-52 Hz frequency range over which the HBO 
was detected was divided in 16 selection bins with 
widths of 2 or 4 Hz, depending on the signal to noise 
level. For each selection interval the power spectra 
were averaged, and a mean S z value was determined. 
The HBO selection criteria were applied to all data 
along the HB and near the hard vertex. When the 
source was near the hard vertex, on the NB, or the 
FB, we selected the power spectra according to the S z 
value. An overlap between the two methods occurred 
for data near the hard vertex; both selection meth- 
ods yielded the same results for the fit parameters to 
well within the statistical errors (see Section |3|). Sep- 
arately, for each set of observations (A,B, and C) we 
also determined the kHz QPO properties according 
to the S z method. 

3. RESULTS 

Using the fit function described by Jonker et al. 
(1998) which consisted of two Lorentzians to describe 
the HBO and the second harmonic of the HBO, and 
a cut-off power law to describe the LFN noise compo- 
nent, we obtained poor fits. Compared with Jonker 
et al. (1998) we combined more data, resulting in 
a higher signal to noise ratio. First we included 
a peaked noise component (called sub-HBO compo- 
nent) at frequencies below the HBO, since a similar 
component was found by van der Klis et al. (1997) 
in Sco X-l. This improved the Xred °f the Re- 
maining problems were that the frequency of the sec- 
ond harmonic was not equal to twice the HBO fre- 
quency (similar problems fitting the power spectra 
on the HB of Cyg X-2 were reported by Kuulkers, 
Wijnands, & van der Klis 1999), and the frequency 
of the sub-HBO component varied erratically along 
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Fig. 1. — Hardness-intensity diagrams for observations 1 and 9-18 (A), 2-8 (B), and 19-24 (C) (see Table 0). The hard color 
is defined as the logarithm of the 9.7-16.0/6.4-9.7 keV count rate ratio. The intensity is defined as the three-detector count rate 
measured in the 2-16.0 keV band. The data were background subtracted but no dead-time correction was applied. The dead-time 
correction factor was less than 1.5%. 
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Fig. 2. — The S z value of ~150 individual 16 s length power spectra from observation 1 plotted against their fitted HBO 
frequency. The line represents the best linear fit. The Xred is 1-75 for 144 degrees of freedom, the linear-correlation coefficient is 
0.83. 
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5000 1(T 

Time (s) 

Fig. 3. — The dynamical power spectrum of part of observation 1 showing the 1/16-60 Hz range, with a frequency resolution of 
0.5 Hz in the energy band 5-60 keV. The grey scale represents the Leahy normalized power (Leahy et al. 1983). Data gaps have 
been omitted for clarity. Clearly visible is the HBO at ~ 18 to ~ 50 Hz and the LFN component at low (< 10 Hz) frequencies. 
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the HB. Inspecting the fit showed that both the fit 
to the high frequency tail of the HBO, and the fit to 
its second harmonic did not represent the data points 
very well. Including an additional component in the 
fit function representing the high frequency tail of the 
HBO (called shoulder component after Belloni et al. 
[1997] who used this name) resulted in a better fit 
to the HBO peak, a centroid frequency of the HBO 
second harmonic more nearly equal to twice the HBO 
frequency, and a more consistent behavior of the fre- 
quency of the sub-HBO component (which sometimes 
apparently fitted the shoulder when no shoulder com- 
ponent was present in the fit function). 
We also experimented with several other fit function 
components to describe the average power spectra 
which were used by other authors to describe the 
power spectra of other LMXBs. Using a fit func- 
tion built up out of a broken power law, to fit the 
LFN component, and several Lorentzians to fit the 
QPOs after Wijnands et al. (1999a) results in signif- 
icantly higher Xred va l ues than when the fit function 
described in Section |2| was used (xled= L66 for 205 
degrees of freedom (d.o.f.) versus a Xred = 1-28 with 
204 d.o.f.). We also fitted the power spectra using the 
same fit function as described in Section ||| but with 
the frequency of the sub-HBO component fixed at 
Hz, in order to test whether or not an extra LFN-like 
component centred around Hz was a good represen- 
tation of the extra sub-HBO component. Finally, we 
tested a fit function built up out of two cut-off power 
laws; one describing the LFN and one either describ- 
ing the sub-HBO component or the shoulder compo- 
nent, and three Lorentzians, describing the HBO, its 
second harmonic, and either the sub-HBO or shoul- 
der component when not fitted with the cut-off power 
law. But in all cases the Xred vames obtained using 
these fit functions were significantly higher (for the 
24-26 Hz selection range values of 1.52 for 205 d.o.f., 
1.62 for 204 d.o.f., and 2.00 for 205 d.o.f. were ob- 
tained, respectively). 

Settling on the fit function already described in Sec- 
tion ^j, we applied an F-test (Bevington &: Robinson 
1992) to the x 2 °f the fits with and without the extra 
Lorentzian components to test their significance. We 
derived a significance of more than 8 a for the sub- 
HBO component, and a significance of more than 6.5 
a for the shoulder component, in the average selected 
power spectrum corresponding to HBO frequencies of 
24 to 26 Hz. In Figure | we show the contribution of 
all the components used to obtain the best fit in this 
power spectrum. 

The properties of all the components used in describ- 
ing the low-frequency part of the average power spec- 
tra along the HB are given in Fig. |as a function of 
S z . When the HBO frequency was higher than 32 



Hz, the sub-HBO and shoulder component were not 
significant. We therefore decided to exclude these 
two components from the fit function in the HBO 
frequency selections of 32 Hz and higher. When this 
affected the parameters determined for the remaining 
components in the fit function, we mention so. Split- 
ting the total counts into different photon energies 
reduced the signal to noise in each energy band and 
therefore these effects were more important in the fits 
performed to determine the energy dependence of the 
parameters. 

3.1. The LFN component 

The fractional rms amplitude of the LFN decreased 
function of S z (Fig |5| A), with values ranging 
from 10% to 2.2% (5.0-60 keV). Upper limits on the 
LFN component were calculated by fixing the power 
law index at 0.0. The power law index of the LFN 
component increased from ~ at S z ~ 0.5 to ~ 0.4 
around S z = 0.9; when the source moved on to the 
NB the index of the power law decreased to values 
slightly below 0.0 (Fig || B). The cut-off frequency 
of the LFN component increased clS ct function of S 7 . 
For S z > 1.0 the cut-off frequency could not be de- 
termined with high accuracy (Fig |H| C). 
The LFN fractional rms amplitude depended strongly 
on photon energy all across the selected frequency 
range. The rms amplitude increased from 5% at 2- 
5.0 keV to more than 15% at 13.0-60 keV (S z =0.48). 
The power law index, a, of the LFN component was 
higher at lower photon energies (changing from 0.3- 
0.5 along the HB at 2-5.0 keV) than at higher pho- 
ton energies (changing from —0.2-0.2 along the HB 
at 13.0-60 keV). The cut-off frequency of the LFN 
component did not change as a function of photon 
energy. 

3.2. The HBO component 

The fractional rms amplitude of the HBO decreased 
function of S z (Fig || G), with values ranging 
from 10% to 1.7% over the detected range (5.0-60 
keV). Upper limits on the HBO component were de- 
termined using a fixed FWHM of 15 Hz. The fre- 
quency of the HBO increased function of S 7 
but for S z >1.0 it was consistent with being constant 
around 50 Hz (Fig|l). 

The ratio of the rms amplitudes of the LFN and the 
HBO component, of interest in beat frequency mod- 
els (see Shibazaki & Lamb 1987) decreased from ~ 1 
at an S z value of 0.48 to ~ 0.6 at S z values of 0.8- 
1.0. The ratio increased again to a value of ~ 0.9 at 
S z =1.05 when the source was on the NB. 
The HBO rms amplitude depended strongly on pho- 
ton energy all across the selected frequency range. 
The rms amplitude increased from 5% at 2-5.0 keV 
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Fig. 4. — Leahy normalized power spectrum showing the different components used to fit the 5.0-60 keV power spectrum. 
The full line represents the LFN component, and the constant arising in the power spectrum due to the Poisson counting noise; 
the dashed line represents the sub-HBO component; the dotted line represents the HBO; the dashed-dotted line represents the 
shoulder component; and the dash-three dots-dash line represents the harmonic of the HBO. 
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Fig. 5. — (A) Rms amplitude of the low-frequency noise (LFN); (B) power law index of the LFN; (C) cut-off frequency of the 
LFN; (D) rms amplitude of the noise component at frequencies below the HBO frequency (sub-HBO component); (E) FWHM of 
the sub-HBO component; (F) frequency of the sub-HBO component; (G) rms amplitude of the HBO; (H) FWHM of the HBO; 
(I) frequency of the HBO; (J) rms amplitude of the shoulder component used to describe the HBO; (K) FWHM of the shoulder 
component; (L) frequency of the shoulder component; (M) rms amplitude of the harmonic of the HBO; (N) FWHM of the har- 
monic; (O) frequency of the harmonic. The points represent data selected according to the HBO selection method and the bullets 
represent the data selected according to the S z selection method (parameters measured in the 5.0-60 keV band, see text). The 
two methods overlap starting around S z ~1.0. 
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to ~ 16% at 13.0-60 keV (at S z =0.48) (see Figure g 
[dots] for the HBO energy dependence in the 26- 
28 Hz range). The increase in fractional rms am- 
plitude of the HBO towards higher photon energies 
became less as the frequency of the HBO increased. 
At the highest HBO frequencies the HBO is relatively 
stronger in the 8.4-13.0 keV band than in the 13.0-60 
keV band. The ratio between the fractional rms am- 
plitude as a function of photon energy of the HBO at 
lower frequencies and the fractional rms amplitude as 
a function of photon energy of the HBO at higher fre- 
quencies is consistent with a straight line with a pos- 
itive slope. The exact fit parameters depend on the 
HBO frequencies at which the ratios were taken. This 
behavior was also present in absolute rms amplitude 
(= fractional rms amplitude *I X , where I x is the count 
rate), see Fig. |7| So, this behavior is caused by actual 
changes in the QPO spectrum, not by changes in the 
time-averaged spectrum by which the QPO spectrum 
is divided to calculate the fractional rms spectrum of 
the HBO. The FWHM and the frequency of the HBO 
were the same in each energy band. 

3.3. The second harmonic of the HBO 

The rms amplitude of the second harmonic of the 
HBO decreased as a function of S z (Fig [5] M) from 
5.2% to 3.6% (5.0-60 keV). Upper limits on the sec- 
ond harmonic of the HBO were derived using a fixed 
FWHM of 25 Hz. The frequency of the second har- 
monic of the HBO was consistent with being twice the 
HBO frequency when the sub-HBO and the shoulder 
component were strong enough to be measured (see 
Fig H O, and Fig. ||). When these two extra compo- 
nents could not be determined significantly, due to 
the limited signal to noise, and we therefore omit- 
ted them from the fit function (as explained above), 
the frequency of the second harmonic of the HBO 
was clearly less than twice the HBO frequency (see 
Fig. |). 

The rms amplitude of the second harmonic of the 
HBO was also energy dependent. Its rms amplitude 
increased from less than 4% in the 2-5.0 keV, to more 
than 9% in the 8.6-13 keV band. The FWHM of the 
second harmonic varied erratically in the range of 10— 
50 Hz. This is not necessarily a property of the second 
harmonic since the HBO shoulder component which 
was not significant by themselve was omitted from 
the fit function. This may have influenced the fit to 
the FWHM of the second harmonic when it was weak. 
Its frequency was consistent with being the same in 
each energy band. 

3.4. The sub-HBO component 

The centroid frequency (Fig [5| F) and FWHM (Fig || 
E) of the Lorentzian at sub-HBO frequencies in- 



creased from 9.3±0.3 Hz to 13.6±1.0 Hz and from 
7.1±0.6 Hz to 18±3 Hz, respectively, as the source 
moved up the HB from S z = 0.48 to 0.67. The rms 
amplitude of this component did not show a clear re- 
lation with S z ; its value was consistent with being 
constant around 5% (Fig [| D). Upper limits on the 
sub-HBO component were determined using a fixed 
FWHM of 15 Hz. The frequency of the sub-HBO 
component is close to half the frequency of the HBO 
component. The fact that the ratio between the HBO 
frequency and the sub-HBO frequency is not exactly 
2 but ~ 2.2 may be accounted for by the complexity 
of the data and therefore its description. 
We detected the sub-HBO component in the three 
highest energy bands that we defined over an S z range 
from 0.48-0.65. Its rms amplitude is higher in the 
highest energy band (~7% in 13.0-60 keV, and less 
than 5% in 6.4-8.6 keV) and decreased as a function 
of S z , while the FWHM and the frequency increased 
from 6-12 Hz, and 9-15 Hz, respectively. 



3.5. The HBO shoulder component 

At an S z value of 0.48 (the left most part of the HB) 
the frequency of the shoulder component was higher 
than the frequency of the HBO, and the frequency 
separation between them was largest (Fig [5] L and 
Fig. |8|). Both the frequency of the shoulder and the 
HBO increased when the source moved along the HB, 
but the frequency difference decreased. The FWHM 
of the shoulder component increased from 7±2 Hz 
to 17±3 Hz and then decreased again to 8.6±1.0 Hz 
as the frequency of the HBO peak increased from 
25.7±0.7 Hz to 30.9±0.4 Hz (Fig § K). From an S z 
value of 0.61 to 0.67 the frequency was consistent with 
being constant at a value of 32 Hz. The rms ampli- 
tude was consistent with being constant around 4% 
(5.0-60 keV), over the total range where this compo- 
nent could be detected (Fig [| J), but the data is also 
consistent with an increase of fractional rms ampli- 
tude with increasing HBO frequency. Upper limits on 
the HBO shoulder component were determined using 
a fixed FWHM of 7 Hz. 

In the various energy bands the HBO shoulder com- 
ponent was detected seven times in total; once in the 
2-5.0 keV band, three times in the 6.4-8.6 keV band, 
and three times in the in the 8.6-13.0 keV band, with 
rms amplitudes increasing from ~3% in the 2-5.0 keV 
band to ~6% in the 8.6-13.0 keV band, and a FWHM 
of ~10 Hz. Upper limits of the order of 3%-4%, and 
of 5%-7% were derived in the two lowest and three 
highest energy bands considered, respectively. These 
are comparable with or higher than the rms ampli- 
tudes of this component determined in the 5-60 keV 
band. 
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Fig. 6. — The figure shows the typical energy dependence of the rms amplitude of the HBO (bullets) and NBO (squares) as 
measured in the frequency range 26-28 Hz for the HBO, and as measured in the S z 1.0-1.9 range for the NBO. 
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a function of the photon energy. 
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Fig. 8. — The frequencies of the four Lorentzian components used to describe the average 5.0-60 keV power spectra, as a 
function of the HBO frequency. Shown are from low frequencies to high frequencies; the sub-HBO component (stars), the HBO 
(bullets), the shoulder component (open circles), and the second harmonic of the HBO (squares). The solid line represents the 
relation v = 0.5 * vhbo, the dashed-dotted line represents v — 1.0 * vhbo, and the dotted line represents v — 2.0 * vhbo- Errors 
in the HBO frequency are in some cases smaller than the symbols. 
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Table 2 

Properties of the NBO fitted using one or two Lorentzians, as a function of S z in the 5.0-60 keV 

BAND. 



"Parameter fixed at this value 



3.6. The NBO component 

The NBOs were not observed when the source was 
on the HB, with an upper limit of 0.5% just before 
the hard vertex (for an S z value of 0.96). They were 
detected along the entire NB and they evolved into 
a broad noise component on the FB. The properties 
arc listed in Table g. The rms amplitude of the NBO 
gradually increased while the source moved from the 
upper NB to the middle part of the NB where the rms 
amplitude is highest. On the lower part of the NB the 
NBO rms amplitude gradually decreased. Upper lim- 
its on the NBO components were determined using a 
fixed FWHM of 5 Hz. 

As the NBO got stronger towards the middle of the 
NB the profile of the NBO became detectably asym- 
metric (see Fig. |9|); between S z =1.25 and 1.54 the 
NBO was fitted using two Lorentzians. The FWHM 
of the NBO as a function of the position along the 
NB was first decreasing from ~10 Hz at S z =1.038 
to around 2.5 Hz on the middle part of the NB (S z 
values from 1.25-1.54), and then increased again to 
~ 5 Hz on the lowest part of the NB. 
Due to the fact that the NBO profiles had to be fit- 
ted using two Lorentzians in part of the data, the 
behavior of the NBO frequency cLS Si function of S 7 
is also not determined unambiguously. Therefore, we 
weighted the frequencies of these two Lorentzians ac- 
cording to one over the square of the FWHM. The 
FWHM weighted average of the two centroid frequen- 
cies of the two Lorentzians used to describe the NBO 
was consistent with a small increase as a function of 
S z from 5.27 ± 0.07 Hz at S 7 =1.25 to 5.83 ± 0.18 Hz 



at S z =1.54. 

We combined all power spectra with an S z between 
1.0 and 1.9 in order to investigate the energy depen- 
dence of the NBO. The rms amplitude of the NBO 
increased as of function of photon energy (see Fig- 
ure |6| [squares] ) . 

3.7. KHz QPOs 

Using the HBO frequency selection method in all data 
combined, the frequency of the kHz QPO peaks in- 
creased from 1971^ Hz to 565_? 4 Hz and from 535±|| 
Hz to 840 ± 21 Hz for the lower and upper peak, re- 
spectively, while the frequency of the HBO increased 
from 20.55 ± 0.02 Hz to 48.15 ± 0.08 Hz. Using the 
S z selection method on the three data sets we defined 
in Section ^ (Figure [[]), we found that the relation 
between the kHz QPO and the HBO is consistent 
with being the same in all three data sets (Fig. 1C 
upper panel). The same relation was found when we 
combined all data and selected the power spectra ac- 
cording to the HBO frequency. 

Upper limits on the kHz QPOs were determined with 
the FWHM fixed at 150 Hz. When only one of the 
two kHz QPO peaks was detected the upper limit 
on the other peak was determined by fixing the fre- 
quency at the frequency of the detected peak plus or 
minus the mean difference frequency between the two 
peaks, depending on whether the lower or the upper 
peak was detected. The properties of the kHz QPOs 
as determined in all data combined when selected ac- 
cording to the HBO frequency are listed in Table |. 
The kHz QPO peak separation was consistent with 
being constant at 339±8 Hz over the observed kHz 
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Fig. 9. — Typical Leahy normalized power spectrum on the NB showing the NBO in the energy band 5.0-60 keV (in the S z 
1.2-1.5 range). The asymmetry of the profile is clearly visible; the drawn line represents the best fit model, using two Lorentzian 
peaks. The dotted line and the dash-dotted line represent the two Lorentzians. 
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Fig. 10. — Upper panel: Relation between the lower and upper kHz QPO peak frequencies and the HBO frequency, as measured 
using all the data selected according to their HBO frequency in the 5-60 keV energy band (filled large squares), and using data 
from observations 1, and 9-18 combined (bullets; see Jonker et al. 1998), data from observations 2-8 combined (stars), and data 
from observations 19-25 combined (diamonds) selected according to the S z selection method. The error bars on the HBO frequency 
are small compared to the size of the data points, and are therefore omitted. Lower panel: The peak separation vs. upper kHz 
QPO frequency as measured when selected according to HBO frequency. The solid, dashed, and dash-dotted lines represent the 
predicted relations between the peak separation and the Keplerian frequency in the Stella & Vietri (1999) model for a neutron 
star mass of 1.4, 2.0, and 2.2 M , respectively. 
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Table 3 

Properties of the kHz QPOs as determined in all 5.0-60 keV data combined, selected according to 

the hbo frequency. 



"Parameter fixed at this value 



QPO range (Fig. [H] lower panel), but a decrease 
towards higher upper peak frequencies similar to 
that found in Sco X-l (van der Klis et al. 1997), 
4U 1608-52 (Mendez et al. 1998), 4U 1735-44 (Ford 
et al. 1998), 4U 1702-429 (Markwardt, Strohmayer, 
& Swank 1999), and 4U 1728-34 (Mendez k van der 
Klis 1999) cannot be excluded. The FWHM of nei- 
ther the lower nor the higher frequency kHz QPO 
peak showed a clear relation with frequency. The rms 
amplitude of the lower and upper kHz QPO peak de- 
creased from 3.1% to 1.1%, and from 4.2% to 1.2%, 
respectively when the HBO frequency increased from 
20.55 to 48.15 Hz. 

4. DISCUSSION 

In the present work we combined all RXTE data 
presently available for the Z source GX 340±0 us- 
ing our new selection method based on the frequency 
of the HBO peak. This allowed us to distinguish 
two new components in the low-frequency part of the 
power spectrum. 

These two extra components were strongest when the 
source was at the lowest count rates on the HB (see 
Fig. |), between S z = 0.48-0.73, i.e., at the lowest 
inferred M. The frequency of one of these compo- 
nents, the sub-HBO component, is close to half the 
frequency of the HBO component. The frequency ra- 
tio was consistent with being constant when the fre- 
quency of the sub-HBO changed from 9 to 14 Hz. A 
similar feature at sub-HBO frequencies has been re- 
ported by van der Klis et al. (1997) in Sco X-l. Since 
the frequency of this component is close to twice the 
predicted Lense-Thirring (LT) precession frequency 
for rapidly rotating neutron stars (Stella & Vietri 



1998), we shall discuss the properties of this com- 
ponent within this framework. 

The other component we discovered, the HBO shoul- 
der component, was used to describe the strong excess 
in power in the HBO profile towards higher frequen- 
cies. If this shoulder component is related to the HBO 
and not to a completely different mechanism which 
by chance results in frequencies close to the frequency 
of the HBO, it can be used to constrain the forma- 
tion models of the HBO peak. We demonstrated that 
both the HBO and the NBO have a similar asymmet- 
ric profile. In the NBO this was previously noted by 
Priedhorsky et al. (1986) in Sco X-l. We shall con- 
sider the hypothesis that the formation of this shoul- 
der is a common feature of the two different QPO 
phenomena, even if the two peaks themselves perhaps 
occur due to completely different physical reasons. 
Our results on the kHz QPOs based on more exten- 
sive data sets at three different epochs and using the 
new HBO selection method are consistent with those 
of Jonker et al. (1998). We discuss the properties 
of the kHz QPOs within the framework of precess- 
ing Keplerian flows (Stella & Vietri 1999), the sonic 
point model (Miller, Lamb, & Psaltis 1998), and 
the transition layer model described by Osherovich 
& Titarchuk (1999), and Titarchuk, Osherovich, & 
Kuznetsov (1999). 

4.1. Comparison with other observations 

In various LMXBs, QPOs have been found whose pro- 
files are clearly not symmetric. Belloni et al. (1997) 
showed that for the black hole candidate (BHC) 
GS 1124-68 the QPO profiles are asymmetric, with 
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a high frequency shoulder. Dieters et al. (1998) re- 
ported that the 2.67 Hz QPO of the BHC 4U 1630-47 
was also asymmetric with a high frequency shoulder. 
In the Z source Sco X-l the NBO profile was also 
found to be asymmetric (Priedhorsky et al. 1986). 
It is clear that asymmetric shapes of the QPO pro- 
files are frequently observed in LMXBs and are not 
restricted to either the black hole candidates or the 
neutron star systems. 

In the BHCs GS 1124-68 (Belloni et al. 1997) and 
XTE J1550-564 (Homan et al. 1999) several QPOs 
were discovered which seem to be harmonically re- 
lated in the same way as we report for GX 340+0, 
i.e. the third harmonic is not detected, while the 
first, the second and the fourth harmonic are. If 
this implies that these QPOs are the same, mod- 
els involving the magnetic field of the neutron star 
for their origin could be ruled out. The time lag 
properties of the harmonic components of the QPOs 
in XTE J 1550-564 are complex and quite distinc- 
tive (Wijnands, Homan, & van der Klis 1999). In 
GX 340+0 no time lags of the harmonic components 
could be measured, but the time lags measured in the 
HBO in the similar Z source GX 5-1 (Vaughan et al. 
1994) are quite different. 

In order to study in more detail the relationship found 
by Wijnands & van der Klis (1999) between the QPOs 
and the noise break frequency in the power spectrum 
of LMXBs, we fitted the LFN component using a 
broken power law. To determine the value for the 
break frequency we fixed the parameters of all other 
components to the values found when using a cut- 
off power law to describe the LFN. Wijnands <fe van 
der Klis (1999a) reported that the Z sources did not 
fall on the relation between the break and QPO fre- 
quency established by atoll sources and black hole 
candidates. They suggested that the Z source LFN 
is not similar to the atoll HFN but the noise com- 
ponent found in Sco X-l at sub-HBO frequencies is. 
By using the centroid frequency of that peaked noise 
component as the break frequency instead of the LFN 
break frequency, the HBO frequencies did fall on the 
reported relation. On the other hand, we find that 
using the sub-HBO frequency instead of the HBO fre- 
quency together with the LFN break frequency, the 
Z source GX 340+0 also falls exactly on the relation. 
Therefore, the suggestion made by Wijnands & van 
der Klis (1999a) that the strong band-limited noise in 
atoll and Z sources have a different origin is only one 
of the two possible solutions to the observed discrep- 
ancy. Our proposed alternative solution is that the 
Z and atoll noise components are the same, but that 
it is the sub-HBO in Z sources which corresponds to 
the QPO in atoll sources. An argument in favour 
of the noise components in Z and atoll sources be- 



ing the same is that the cut-off frequency of the LFN 
component increased as a function of S z , in a simi- 
lar fashion as the frequency associated with the atoll 
high frequency noise (van der Klis 1995, Ford & van 
der Klis 1998, van Straaten et al. 1999). 
Following Psaltis, Belloni, & van der Klis (1999) we 
plotted the sub-HBO frequency against the frequency 
of the lower-frequency kHz QPO. The sub-HBO does 
not fall on the relation found by Psaltis, Belloni, & 
van der Klis (1999) between the frequency of the HBO 
and the lower-frequency kHz QPO frequency. Instead 
the data points fall between the two branches defined 
by the HBO-like QPO frequencies vs. the lower kHz 
QPO frequency at high frequencies (see Psaltis, Bel- 
loni, & van der Klis 1999). 

4.2. HBO - kHz QPO relations 

4.2.1. Lense-Thirring precession frequency 

Stella & Vietri (1998) recently considered the possi- 
bility that the HBO is formed due to the LT preces- 
sion of the nodal points of sligthly tilted orbits in the 
inner accretion disk, but as they already mentioned 
the Z sources GX 5-1 and GX 17+2 did not seem to 
fit in this scheme. For reasonable values of I/M, the 
neutron star moment of inertia divided by its mass, 
the observed frequencies were larger by a factor of 
~2 than the predicted ones. Jonker et al. (1998) 
showed that for GX 340+0 the predicted frequency 
is too small by a factor of 3, if one assumes that the 
higher frequency peak of the kHz QPOs reflects the 
Kcplerian frequency of matter in orbit around the 
neutron star, and that the mean peak separation re- 
flects the neutron star spin frequency. Using the same 
assumptions Psaltis et al. (1999) also concluded that 
a simple LT precession frequency model is unable to 
explain the formation of HBOs in Z sources. 
Detailed calculations of Morsink & Stella (1999) even 
worsen the situation, since their calculations lower 
the predicted LT frequencies. They find that the 
LT precession frequencies are approximately a fac- 
tor of two too low to explain the noise compo- 
nents at frequencies ~20-35 Hz observed in atoll 
sources (4U 1735-44, Wijnands & van der Klis 1998c; 
4U 1728-34, Strohmayer et al 1996, Ford & van der 
Klis 1998). Stella & Vietri (1998) already put forward 
the suggestion that a modulation can be produced at 
twice the LT precession frequency if the modulation 
is produced by the two points where the inclined or- 
bit intersects the disk plane (although they initially 
used this for explaining the discrepancy of a factor 
of two between the predicted and the observed LT 
precession frequencies for the Z sources). 
The sub-HBO peaked noise component we discovered 
could be harmonically related to the HBO compo- 
nent. If the sub-HBO is the second harmonic of the 
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fundamental LT precession frequency, as needed to 
explain the frequencies in the framework of the LT 
precession model where the neutron star spin fre- 
quency is approximately equal to the frequency of 
the kHz QPO peak separation, the HBO must be the 
fourth and the harmonic of the HBO must be the 
eighth harmonic component, whereas the sixth and 
uneven harmonics must be much weaker. This poses 
strong (geometrical) constraints on the LT precession 
process. On the other hand, if the HBO frequency is 
twice the LT precession frequency, which implies a 
neutron star spin frequency of ~900 Hz (see Morsink 
& Stella 1999), the frequency of the sub-HBO com- 
ponent is the LT precession frequency, and the fre- 
quency of the second harmonic of the HBO is four 
times the LT precession frequency. In that case only 
even harmonics and the LT precession frequency are 
observed. 

4.2.2. Magneto spheric beat frequency and 
radial-flow models 

In this section, we discuss our findings concerning 
the QPOs and the LFN component in terms of the 
magnetic beat frequency model where the QPOs are 
described by harmonic series (e.g. Shibazaki & Lamb 
1987). 

If the sub-HBO frequency is proven not to be har- 
monically related to the HBO, the sub-HBO peak 
might be explained as an effect of fluctuations enter- 
ing the magnetospheric boundary layer periodically. 
Such an effect will be strongest at low HBO frequen- 
cies since its power density will be proportional to 
the power density of the LFN (Shibazaki & Lamb 
1987). If it is the fundamental frequency and the 
HBO its first overtone then the magnetospheric beat 
frequency model proposed to explain the HBO for- 
mation (Alpar & Shaham 1985; Lamb et al. 1985) is 
not strongly constrained. 

Within the beat frequency model the high frequency 
shoulder of the HBO peak can be explained as a sign 
of radial drift of the blobs as they spiral in after cross- 
ing the magnetospheric boundary layer (Shibazaki 
& Lamb 1987). Shibazaki & Lamb (1987) describe 
another mechanism which may produce a high fre- 
quency shoulder. Interference between the LFN and 
the QPO caused by a non uniform phase distribution 
of the blobs will also cause the QPO to become asym- 
metric. This effect will be strongest when the LFN 
and the QPO components overlap, as observed. Fi- 
nally, an asymmetric initial distribution of frequen- 
cies of the blobs when entering the magnetospheric 
boundary layer may also form an asymmetric HBO 
peak. 

The changes in the power law index of the LFN as a 
function of photon energy can be explained by vary- 



ing the width or the steepness of the lifetime distri- 
bution of the blobs entering the magnetic boundary 
layer (Shibazaki & Lamb 1987). The decrease in in- 
crease of both the fractional and absolute rms am- 
plitude of the HBO clS Si function of energy towards 
higher frequencies (Fig. 0) also constrains the detailed 
physical interactions occurring in the boundary layer. 
Fortner et al. (1989) proposed that the NBO is 
caused by oscillations in the electron scattering op- 
tical depth at the critical Eddington mass accretion 
rate. How a high frequency shoulder can be produced 
within this model is not clear. Both the HBO and 
the NBO shoulder components were detected when 
the rms amplitude of the HBO and the NBO was 
highest. In case of the NBO, this may be a result 
of the higher signal to noise. Since the rms ampli- 
tude of the NBO shoulder component is consistent 
with being ~2/3 of the NBO rms amplitude (see Ta- 
ble |^) , combining more observations should increase 
the range over which this shoulder component is de- 
tected, if this ratio is constant along S z . In case of the 
HBO the two components seem to merge. While the 
fractional rms amplitude of the HBO shoulder com- 
ponent increased that of the HBO decreased. When 
the fractional rms amplitudes were comparable, the 
HBO was fitted with one Lorentzian. The rms am- 
plitude of both shoulder components increased in a 
similar way as the rms amplitudes of the NBO and 
the HBO with photon energy. So, the formation of 
these shoulder components seems a common feature 
of both QPO forming mechanisms. 

4.2.3. Radial oscillations in a viscous layer 

In Sco X— 1, Titarchuk, Osherovich, & Kuznetsov 
(1999) interpreted the extra noise component in the 
power spectra (van der Klis et al. 1997) as due 
to radial oscillations in a viscous boundary layer 
(Titarchuk &: Osherovich 1999). If the noise com- 
ponent in Sco X-l is the sub-HBO component in 
GX 340+0, the model of Titarchuk & Osherovich 
(1999) can be applied to the frequencies and de- 
pendencies we found for the sub-HBO component 
in GX 340+0. Fitting our data to the relation be- 
tween the frequency of the extra noise component and 
the Keplerian frequency, using the parameters and 
parametrization given by Titarchuk, Osherovich, & 
Kuznetsov (1999), we obtained a value of CV = 15 for 
GX 340+0. This value is much larger than the value 
obtained for Sco X-l (9.76). According to Titarchuk 
& Osherovich (1999) a higher Cm value implies a 
higher viscosity for the same Reynold's number. 

4.3. KHz QPOs and their peak separation 

Recently, Stella & Vietri (1999) have put forward a 
model in which the formation of the lower kHz QPO 
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is due to the relativistic periastron precession (ap- 
sidal motion) of matter in (near) Keplerian orbits. 
The frequency of the upper kHz QPO peaks is the 
Keplerian frequency of this material. The peak sepa- 
ration is then equal to the radial frequency of matter 
in a nearly circular Keplerian orbit, and is predicted 
to decrease as the Keplerian frequency increases and 
approaches the predicted frequency at the marginally 
stable circular orbit. This model can explain the 
decrease in peak separation as observed in various 



sources (see Section 3.7). 

Beat frequency models stating that the upper kHz 
QPO peak is formed by Keplerian motion at a pre- 
ferred radius in the disk (e.g. the sonic point radius, 
Miller, Lamb, & Psaltis 1998), whereas the lower kHz 
QPO peak formed at the frequency of the beat be- 
tween the neutron star spin frequency and this Ke- 
plerian frequency, cannot in their original form ex- 
plain the decrease in peak separation in these two 
sources. A relatively small extension of the model 
(Lamb, Miller, & Psaltis 1998) can, however, produce 
the observed decrease in peak separation. 
Osherovich & Titarchuk (1999) developed a model in 
which the kHz QPOs arise due to radial oscillations 
of blobs of accreting material at the magnetospheric 
boundary. The lower kHz QPO frequency is in their 
model identified with the Keplerian frequency. Be- 
sides this QPO two eigenmodes are identified whose 
frequencies coincide with the upper kHz QPO peak 
frequency and the frequency of the HBO component 
in the power spectra of Sco X-l (Titarchuk & Os- 
herovich 1999). Interpreting our findings within this 
framework did not result in stringent constraints on 



the model. 

We found that the peak separation is consistent with 
being constant (Fig. [H] A and B), but neither a de- 
crease towards higher M as in Sco X-l, 4U 1608- 
52, 4U 1735-44, 4U 1702-429, and 4U 1728-34 nor 
a decrease towards lower M, as predicted by Stella 
& Vietri (1999) can be ruled out. If the model of 
Stella & Vietri turns out to be the right one the 
mass of the neutron star most likely is in the range 
of 1.8 to 2.2 M (see Fig. [10| B). This is in agree- 
ment with the mass of Cyg X-2 derived by Orosz & 
Kuulkers (1999), and with the masses of the neutron 
stars derived when interpreting the highest observed 
kHz QPO frequencies as due to motion at or near the 
marginally stable orbit (Kaaret, Ford, & Chen 1997; 
Zhang, Strohmayer, & Swank 1997). 
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